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The goal ol quantum benchmarking is to certify that imperfect quantum communication devices 
(e.g., quantum channels, quantum memories, quantum key distribution systems) can still be used 
for meaningful quantum communication. However, the test states used in quantum benchmarking 
experiments may be imperfect as well. Many quantum benchmarks are only valid for states which 
match some ideal form, such as pure states or Gaussian states. We outline how to perform quantum 
benchmarking using arbitrary states of light. We demonstrate these results using real data taken 
from a continuous- variable quantum memory. 



I. INTRODUCTION 

Quantum communication is poised to be the first real- 
world application of quantum information theory. Im- 
pressive efforts are underway to build advanced quantum 
devices, such as quantum channels [1], quantum memo- 
ries [2], and quantum key distribution systems [3], which 
will make long-distance quantum communication a real- 
ity. But our control over the quantum world is still im- 
perfect; it is difficult to create devices whose operation 
faithfully matches the ideal. Internal flaws and exter- 
nal influences can negatively impact a device's operation. 
Since quantum effects are very fragile, there is little toler- 
ance for such deficiencies. Fortunately, it is still possible 
for imperfect devices to facilitate quantum communica- 
tion. 

The task of quantum benchmarking is to certify a de- 
vice to be within the "quantum domain," the operational 
regime where an imperfect device is still useful for quan- 
tum communication. In practice, this is done by prepar- 
ing some test states, probing the device in question with 
these states, and making measurements on the result- 
ing output states. The corresponding theoretical task is 
to calculate the threshold (according to some figure(s) 
of merit) between a quantum device and the best com- 
parible classical device. If we can conclude, based on 
the available information, that the device's performance 
beats the benchmark value, then the device is in the 
quantum domain. Of course, there is some interplay be- 
tween the two sides, especially in the choice of test states. 

A number of benchmarks have been developed in re- 
cent years, especially for continuous variable (CV) en- 
codings of light [4-21]. Many of these benchmarks re- 
quire test states which are pure and/or Gaussian. Such 
states are much more amenable to theoretical analysis 
than more general states of light. Now, if we recognize 
that our quantum communication devices are imperfect, 
we must also recognize that our state preparations may 
be imperfect as well (e.g., due to excess noise or other 
flaws in the preparation procedure). This leads to dis- 



crepencies between the idealized benchmarking situation 
and the experimental realization. Using state tomogra- 
phy or some other process to obtain numerical descrip- 
tions of the prepared states, we may find that these states 
are not consistent with the ideal form (e.g., not pure or 
not Gaussian). In this case, is it valid to apply bench- 
marks which assume the ideal form? Furthermore, many 
of the above benchmarks require an infinite ensemble of 
test states. What conclusions can be made when using 
a finite ensemble (as we must always do in practice) if 
the benchmark requires an infinite ensemble? Addressing 
these practical concerns is an important part of quantum 
benchmarking. 

In this work, we present a general benchmarking frame- 
work which is easily adaptable to different testing con- 
ditions. In particular, our framework allows for quan- 
tum benchmarking using arbitrary states, not just states 
with a special form. Previous progress has been made 
toward this goal using entanglement-based benchmarks, 
but these schemes are still limited to special cases: either 
the test ensemble contains only pure states [15-17, 19], or 
it consists of only two mixed states [18] (in either case, the 
states are not required to be Gaussian) . We extend these 
results to any finite number of mixed states, allowing us 
to derive quantum benchmarks for arbitrary finite ensem- 
bles of test states. There is some freedom allowed by this 
extension, and we outline how to strengthen benchmarks 
to give the largest quantum domain, i.e., the best chance 
of certifying imperfect quantum devices. We then ap- 
ply our benchmarking criterion to an actual continuous- 
variable memory. 

The remainder of this paper is laid out as follows. In 
Sec. II we outline the necessary theory. We show how 
to generalize previous results to an arbitrary number of 
test states and how to optimize this generalization. As an 
application, we consider phase-symmetric test ensembles, 
and derive a simple standard form for the bipartite ma- 
trices used to represent such ensembles in entanglement- 
based benchmarking. In Sec. Ill, we test our general 
benchmarking method using data from a real continuous- 
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variable memory. Finally, in Sec. IV, we make some 
concluding observations. 



II. ENTANGLEMENT-BASED QUANTUM 
BENCHMARKING 

A. Benchmarking background 

We picture any potential quantum communication de- 
vice as a channel which takes quantum states as input 
and gives quantum states as output. The distinction be- 
tween quantum and classical communication devices is 
that classical devices are equivalent to measure and pre- 
pare (MP) channels. In these channels, the input state is 
measured, the measurement outcome is communicated or 
stored in a solely classical manner, and an output state 
is prepared based on this classical data. Any device that 
is equivalent to a MP channel is not a true quantum 
communication device, since communication through MP 
channels is equivalent to classical communication. It is 
often convenient to picture benchmarking in an adversar- 
ial way. In this scenario, we allow the possibility that a 
quantum device could have been replaced by a MP chan- 
nel trying to pass our quantum benchmarking tests. Only 
by conclusively ruling out MP operation can we certify 
a quantum device as genuine. Devices which can be dis- 
criminated from classical channels are said to be in the 
quantum domain. 

In practice, we probe the device with a finite number 
of different test states {pk}^^ 1 and make measurements 
on the resulting output states. From the measurement 
results and knowledge of the test ensemble, we must de- 
termine whether the device is in the quantum domain 
or whether the same operation could be performed by a 
classical channel. With enough test states and measure- 
ments, we can characterize a device exactly, i.e., perform 
process tomography, allowing a device's performance to 
be tested against any theoretical benchmark [22, 23]. 
However, process tomography can be resource intensive, 
especially for CV systems. Our goal will be to certify 
quantum devices with minimal requirements. If a given 
benchmarking test fails, it may still be possible to cer- 
tify a quantum device by incorporating more test states 
and/or more measurements. 

There are a number of approaches for determining 
whether a device is in the quantum domain. Typi- 
cally, some specific operational or theoretical quantities 
are chosen and the corresponding quantum domain is 
mapped out by determining the limits of classical devices. 
One such example is the state-independent T-V bench- 
mark based on signal-transfer and conditional-variance 
[12, 24], parameters which capture how well a signal is 
preserved and how much extra noise is introduced by the 
device. Another much-studied figure of merit is the aver- 
age fidelity between the test states and the corresponding 
output states after the device action (represented by the 



map A): 

M-l 

F=Y / PkF(p k ,A[ Pk }), (1) 

fc=0 

where pk are some fixed weights. If the test states are not 
perfectly distinguishable, a MP channel cannot achieve 
perfect fidelity. Indeed, for a given test ensemble, there 
is some maximum average fidelity F MP which is achiev- 
able by MP channels. A device must beat this benchmark 
value, F > F MP , to be in the quantum domain. This 
approach requires us to compute the optimal value t MP 
and to calculate or bound the experimentally-achieved 
value F. Typically, finding the experimental fidelity re- 
quires tomography on the output states. As well, tomog- 
raphy on the prepared states is needed to ensure that 
they match the required form. 

A number of fidelity-based benchmarks have been 
derived for CV test states, including coherent states, 
squeezed states, and thermalized states [4, 5, 7-11, 13]. 
Other benchmarks use the same approach but slightly 
different figures of merit [11, 14] or models for classi- 
cal channels [6] . Many of these benchmarks require test 
ensembles with an infinite number of states, although fi- 
delity benchmarks for any two pure states have also been 
calculated [8, 26]. The classical threshold F MP is of- 
ten very difficult to determine, making the fidelity-based 
approach hard to generalize to arbitrary test ensembles. 
This limits the applicability of this method to those en- 
sembles where a benchmark value has been found, and 
we must verify that the experimentally-prepared states 
match the ideal form. As well, without (costly) tomogra- 
phy on the output states, or extra assumptions, it can be 
difficult to properly determine the experimental fidelity 
value F. 

There is another approach to quantum benchmarking 
[15-21] which avoids many of these problems. This ap- 
proach is based on the following important property: MP 
channels always break entanglement. In fact, the set of 
MP channels is equivalent to the set of entanglement- 
breaking channels [27]. By showing that a device pre- 
serves entanglement, we can exclude it from the MP class. 
To accomplish this goal, a bipartite entangled state, built 
from the test states, is considered. For the moment, we 
will focus on pure test states pk = IV'fcXV'fcl) and we will 
generalize to mixed states in the next subsection. We 
imagine the test states are coupled with an ancillary sys- 
tem, giving the following entangled state 

M-l 

l* cnt U-7nEl fc MV^. (2) 

k — 

Here, Alice's states lie in an M— dimensional Hilbert 
space for which {[/c)^}^!^ 1 forms an orthonormal basis. 
We let the device in question act on the A' subsystem of 
this entangled state, mapping it to some output system 
B, while subsystem A is kept isolated. The reduced state 
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Pa remains the same, independent of the channel: 



B. Entanglement-based picture for arbitrary states 



\p A \i) = Ji mi>k) 



(3) 



Aside from a fixed prefactor, the density matrix pa is 
the Gram matrix (matrix of overlaps) for the pure test 
states {IV'fc)}- This structure will be important when we 
generalize to arbitrary states in the next subsection. 

Performing the local projective measurements {A^ := 
I^X^LKlo 1 011 l* cnt )AA' effectively prepares the test 
states {\il>k) A'}^^ 1 a t random. After passing through 
the channel, the output states {p% mt ■= A(|V>kXV'k|)}fcLo 
(on the system B), are measured with some fixed 
set of measurements {0 J B }. Using only these mea- 
surement results and knowledge of the Gram matrix 
Pa, we aim to determine whether the state p A u g — 
[icU ® A] |* ont )(* cnt | j4yl , is entangled. 

One method to certify entanglement is with a witness- 
ing procedure [15-17, 19]. On the other hand, it can be 
useful not only to qualify entanglement, but also to quan- 
tify it. An alternative approach to entanglement-based 
benchmarking [20, 21] is therefore to bound the entangle- 
ment of the final state away from zero using some suit- 
able entanglement measure £ , i.e., to show £(pa u b) > 0. 
One advantage of the quantitative method is that it can 
be used to compare different devices or situations which 
fall in the quantum domain. When we demonstrate our 
benchmarking method in Sec. Ill, we will take the quan- 
titative approach, using tools from [21]. 

The measurements on system A are restricted to the 
projections {n^j-f^ 1 . Performing more general mea- 
surements on system A would prepare superpositions of 
the test states. But these superposition states can sim- 
ply be included in the test ensemble, offering the same 
formalism with a slightly different entangled state for Eq. 
(2). This source replacement scheme allows us to work 
with the test states in practice and consider the entangled 
state l^)^/ as a virtual construct. The Gram matrix pa 
provides the theoretical link between the test states and 
the virtual entangled state. 

Finally, the Gram matrix compactly contains all in- 
formation about the pure test states which is needed for 
the given benchmarking procedure. A full description of 
the test ensemble is not strictly necessary. Of course, 
we require some method for determining the overlap in- 
formation. State tomography is one way to obtain this, 
but there may be other less costly ways to constrain the 
Gram matrix. Similarly, we do not require tomography 
on the output states for entanglement-based benchmark- 
ing. Indeed, previous papers have demonstrated that 
useful quantum benchmarks can be found by measur- 
ing only two conjugate quadratures [1, 15-21, 28]. Thus, 
entanglement-based benchmarking can be very practical 
from an experimental point of view. 



How can the entangled state in Eq. (2) be extended to 
the case where the test states {pk^^o are mixed? Some 
progress was made in this direction in [18], restricting to 
the case of only two test states. Among several potential 
ways to generalize Eq. (2), it was found that using pu- 
rifications of the test states was the most useful. In this 
framework, we imagine that instead of preparing the test 



states, we prepare purifications {iTfc/^/^ 
A" is some purifying system, i.e., 

Tr A "(\Tk){T k \ A , A ,,) = Pk- 



M-l 
k=0 J 



where 



(4) 



From now on, \Tj~) will always refer to a purification of 
Pk- 

To clarify, the experimental procedure (using the test 
states {pk}) remains the same. Theoretically, the bench- 
marking procedure is analyzed as if the purifications 
{IF/Oa'A") were use d instead, i.e., the device is allowed 
to act on the larger system A'A" . Measurements are still 
restricted to the B subsystem, and the only information 
we retain about the purifications for the final benchmark- 
ing are the overlaps (r^ | T;) . In this way, we give extra 
power to an adversarial MP device, since the purifications 
may be more easily distinguished than the test states. Al- 
though this modification may result in a slightly smaller 
quantum domain, any benchmarks based on the purifi- 
cations are valid. Now, even though we give extra power 
to an adversarial device, we should still aim to make 
our quantum benchmarks as stringent as we can, i.e., 
to keep the quantum domain as large as possible. Ref. 
[18] offers one method to limit the adversarial advantage 
provided by purifications, but this introduces additional 
numerical complexity, and we will not consider it in the 
present paper. However, there still remains some free- 
dom in the actual choice of purifications, which we can 
use to strengthen our benchmarking scheme. 

For two mixed test states po an d Pi , the virtual entan- 
gled state takes the form 



I*) 



AA'A' 



7^(l°>Al r oW 



A'A'- 



(5) 



Purifications arc not unique, but we must fix some choice 
to proceed further. Since we only retain information 
about the overlap (Fo|ri), it is important to know the 
range of allowed values for this parameter. We can al- 
ways find purifications which are orthogonal. On the 
other hand, from Uhlmann's theorem [29, 30], we have 



max | (r | Ti> 



Tr 



PoPiVPo 



F(Po,Pi) 



(6) 

where F(p , p{) is the fidelity between the two test states. 
Thus, the range of allowed overlaps is | (r |ri) | 2 £ 
[0, F(p , pi)}. In [18], it was argued that the best bench- 
marks (i.e., those which give the best opportunity to 
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certify devices in the quantum domain) come from pu- 
rifications which saturate the Uhlmann fidelity bound. 
There was numerical support for this claim, but one can 
also appeal to the fact that the fidelity is a measure 
of (in)distinguishability. Accordingly, choosing purifica- 
tions which are as indistinguishable as possible will make 
it harder for a MP channel to simulate a quantum chan- 
nel reliably. We will make this intuition more precise 
in the following pages. Finally, we note that Uhlmann's 
theorem has also been suggested for fidelity-based bench- 
marking using two mixed states [8]. 

It is straightforward to generalize the entangled state 
in Eq. (5) to accomodate M > 2 mixed states. Specifi- 
cally, we will use 

A/-1 

l*W» = ^El fc )Al r fe)A'A»- (7) 
k— 

The reduced density matrix pa — Tta'A" \^)(^\aa'A" * s 
now (aside from the prefactor) the Gram matrix of the 
purifications {iTfe)}*^ 1 : 

(fc|PAlO = ^(Ti|r fc ). (8) 

It is not straightforward how to choose purifications in 
this case. The fidelity is a useful measure for two states, 
but there is no clear extension of the concept for mul- 
tiple states. Indeed, the set of allowed values for the 
overlaps (r fe |r ; ) has a more subtle structure than just a 
product of the ranges allowed by Eq. (6). It may not 
be possible to choose purifications which achieve the fi- 
delity bound for all pairs. In order to determine the best 
choice of purifications, we will examine the structure of 
entanglement-based benchmarking at a deeper level. 



C. Optimizing the purifications 

How can we choose purifications to optimize the 
strength of entanglement-based benchmarks based on 
multiple mixed test states? To begin answering this, we 
make the following observation. 

Proposition 1. Let {Ir*,.)}*^ 1 and {lA^}^ 1 be 
two different sets of purifications of the test states 
{pk\h=o- Assume that there is a completely positive, 
trace-preserving (CPTP) map which transforms one set 
of purifications to the other, i.e., there exists fl such that 

n[|r*xr*|] = |A fc XA fc | (9) 

for all k. If we cannot conclude that a device is in the 
quantum domain with information based on the {|Afc)}, 
then we cannot conclude that the same device is in the 
quantum domain using information about the {iTfe)}- 
Proof. Represent the device under investigation by the 
map A. This map takes in states on the joint A' A" sys- 
tem and outputs states on the B system. We proceed in 



two steps. First, we assume that we can perform tomo- 
graphically complete measurements on the output states. 
Later we will relax this to incomplete measurements. In 
either case, we work under the constraint that, for all k, 

A[\A k ){A k \] = A[\T k ){T k \] = p%»\ (10) 

where p£ ut is the output state corresponding to the test 
state p k . This constraint is experimentally enforced, 
since the output states are independent of the purifica- 
tion we use to theoretically describe the protocol. Con- 
sider the case where we describe the test states using the 
purifications {|Afc)}. A benchmarking protocol is unsuc- 
cessful when we cannot discriminate between the device 
and a MP channel, i.e., when there is a MP channel which 
gives the same output states as the device. Assume then 
that there exists an MP channel Amp such that 

A M p[|A fc )(A fc |] =pT" (11) 

for all k. Now consider the alternate situation where we 
use the purifications {ITfe)}, and there exists a CPTP 
map fl as in Eq. (9). We define another channel £mp 
by 

Smp := A M p ° 0. (12) 

The concatenation of any channel with a MP channel is 
as a MP channel, so £mp is in the MP class. From Eqs. 
(9-11), we must have 

5W[|r*xr fc |] = P r (is) 

for all k. Therefore, there exists a MP channel which 
gives the observed output states when using the purifica- 
tions {|r fc )}. 

Even when we do not have enough measurements for 
complete tomography, we can follow simular arguments. 
Instead of having one output state p£ ut for every k, we 
have a set of states C[p^ ut ], all of which have the same 
expectation values with respect to the employed mea- 
surement operators {Og}: 

C[ P n := {r\Tr(r6i) = (0' B ) mcas V j). (14) 

In this CcLS6, cl benchmarking protocol is unsuccessful 
when there is a MP channel Amp such that 

A M p[|A fc )(A fc |] eC[ P ° fc ut ] (15) 

for all k. Using this channel, we again define a MP chan- 
nel Smp as in Eq. (12). By Eqs. (9) and (15), we 
conclude 

X MP [\T k ){T k \] GC[pH- (16) 

In either case, if benchmarking is unsuccessful for the 
{|Afc)}, it cannot be successful for the {|Ffc)}. □ 
This proposition has some important consequences. 
For one, it tells us that in the source-replacement scheme 
we are using, the quantity of entanglement does not have 
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a direct significance [31]. Instead, the focus should be 
to find those purifications which cannot be collectively 
transformed, via CPTP maps, to any other valid set. 
Benchmarking schemes built with such purifications pro- 
vide the hardest challenge for an adversarial MP device 
attempting to mimic a true quantum device. Indeed, 
such limiting purifications necessarily lead to a larger 
quantum domain than any other comparable choice. Be- 
cause we only use the overlaps for benchmarking, we can 
optimize the CPTP map condition on the level of Gram 
matrices. In the rest of this paper, we will use the terms 
'purifications' and 'Gram matrix (of the purifications)' 
somewhat interchangeably. 

The following theorem will help us translate the CPTP 
map condition to a more amenable form. 

Theorem 1 ([32-34]) Let {\ lh )}™S l and {\5 k )}^ be 
two sets of pure states ( they do not need to be purifica- 
tions). Let G and D be the corresponding Gram matrices, 
with elements Gij := (7j|7i) and Dij := (Sj\Si). There 
exists a CPTP map taking the former states to the 
latter if and only if the Gram matrices are related by 

G = PoD, (17) 

where o denotes the Hadamard (or Schur or entrywise) 
product. The matrix P satisfies P > and its diagonal 
elements are given by diag(P) = {1,1,1,...}. 

To avoid potential confusion, we point out that while 
the CPTP map f2 takes the {|7 fc )} to the {\S k )}, the 
Gram matrix condition, Eq. (17), has the opposite sense, 
i.e., G is obtained by doing a particular operation on D. 
Since some of the matrix elements of P could be zero, we 
cannot invert the equation to give D as a function of G, 
so we leave the relation in this form. We also point out a 
corollary to this theorem: all compatible purifications of 
the test states {pk} can be prepared by applying a CPTP 
map to a set of purifications which are orthogonal. There- 
fore, orthogonal purifications can be seen as generators 
for the rest of the set of purifications. Of course, we are 
interested in the other end of this generation, i.e., the 
limiting sets of purifications. 

Given some ensemble of test states, it may be quite 
difficult to determine the best purifications analytically, 
especially for arbitrary test states. Alternatively, we 
can attempt to find limiting purifications by maximizing 
some objective function /, defined on ensembles of pure 
states (or on the corresponding Gram matrix) , which pre- 
serves the order structure induced by CPTP maps. In 
other words, if there exists a CPTP map £1 taking the 
pure states {iTfc)}^!^ 1 to the pure states {l^fe)};^ 1 , tne 
desired function must satisfy 

/({|7*>}) < /({|4>}). (18) 

When this property holds, the purifications which are 
limiting in the sense of CPTP maps will maximize the 
objective function /. Before discussing candidates for the 



objective function, we point out one caveat. Namely, the 
order structure imposed on the purifications by CPTP 
maps is not a total order. Accordingly, purifications 
which are not linked by a CPTP map are not compara- 
ble. Nevertheless, we can design quantum benchmarking 
schemes with any feasible set of purifications, and this 
heuristic provides a way to obtain strong candidates. 

A good candidate for the objective function is the 
Gram matrix purity V (recall from Eq. (8) that pa is 
essentially the Gram matrix of the chosen purifications) : 

M-l M-l 

fc=o ;=o 

M-l M-l 

£i< r *i r <>i 2 - as) 

k=0 1=0 

From the above theorem, when there is a CPTP map 
taking {Ir^)} — > {\A k )}, we must have 

(r fe |r,) = p kl (A k \A l ), (20) 

with \P k i\ < 1. Therefore the Gram matrix purity is 
monotonic with respect to CPTP maps, as required. 

The Gram matrix purity also has links to the distin- 
guishability of the test states. For example, when work- 
ing with only two test states po and p%, the maximal 
value of V is 

max7>=i(2 + 2 max I (r | Ti) I 2 ) (21) 
4 V |r D >,|r x ) ) 

=i(l + F(p 0) Pi)). (22) 

Aside from a fixed affine transformation, the Gram ma- 
trix purity is the fidelity between the test states. For 
M > 2, the Gram matrix purity defines a kind of aver- 
aged multi-state analog of the fidelity. Finally, consider 
the state 

M-l 

^^M^^U"' ( 23 ) 

fe=0 

Since a device has no information about which test state 
was prepared, this is the effective state input to the device 
in the purification picture. But since the virtual entan- 
gled state in Eq. (7) is pure, the spectra of the reduced 
states are equal, i.e., spec(pA) = spec(pA>A")- Therefore, 
when pa is of high purity, then so is the average test state 
Pa 1 A"- Intuitively, when the average state pa'A" is close 
to a pure state, its constituant test states {|Ffe)(rfe|} will 
be difficult to distinguish. 

One drawback of the Gram matrix purity is that it will 
be difficult to compute the maximum for M > 2, where 
we do not currently have an analytic formula like Eq. (6) . 
Alternatively, since the objective function in Eq. (19) is 
convex and the set of compatible Gram matrices is con- 
vex (see Appendix A), we could maximize this function 
numerically using the methods of convex optimization. 
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Although this presents one path towards our goal, we 
will not pursue it here. 

Instead, the Gram matrix purity helps us motivate a 
slightly different objective function which is easier to op- 
timize numerically. To write this alternate function, we 
abbreviate the overlaps by 

Z u := (T k \Ti) (24) 

and decompose them into real and imaginary parts: 

Z kl =X kl +iY kl . (25) 

The alternate objective function is given by 

M-1 M-1 

h := Y Y ( Xkl + y «) ' ( 26 ) 
k=l+l 1=0 

Aside from a prefactor and duplicated terms, the main 
difference between h and the Gram matrix purity V is 
that the squared modulus of the overlaps is replaced by 
a sum of the real and imaginary parts. The advantage of 
this objective function is that it is linear in the param- 
eters Xki and Y k i, allowing us to numerically optimize 
it using a semidefinite program, which can be done effi- 
ciently In cases where the the imaginary parts Yki van- 
ish, then X k i +Y k i is a lower bound to the modulus \Z k i\. 
As well, when the objective function h is large, the Gram 
matrix purity will also be large, meaning we are close to a 
limiting Gram matrix. The main drawback of this objec- 
tive function is that it is not always monotonic. Never- 
theless, we remind the reader that any compatible Gram 
matrix, including one found by an optimization process 
over h, is valid for benchmarking purposes. Indeed, we 
will see in Sec. Ill that, despite being non-monotonic, 
this function leads to a Gram matrix with near-optimal 
purity. 

Before moving on, we pause to summarize the main 
results of this section. Given some arbitrary ensemble 
of test states {pfc}^.-^ 1 , we consider a virtual entangled 
state as in Eq. (7), where the states {\^k) a'A"}^^ 1 are 
purifications of the test states. Although the test states 
are used in practice, we benchmark a device by assum- 
ing that the device has access to the purifications. This 
framework provides more power to an adversarial device, 
but there is also some flexibility in the choice of purifica- 
tions, which can be used to restrict the adversarial advan- 
tage. We showed that purifications which are extremal 
in the sense of CPTP maps are better for benchmarking 
than all other comparable choices, so we should aim to 
build benchmarks using such purifications. We suggest 
finding good candidate purifications by optimizing some 
appropriate function which is monotonic with respect to 
CPTP maps. We propose the Gram matrix purity V 
as a suitable candidate for this objective function, al- 
though it may be difficult to compute. We also propose 
a numerically simpler, though non-monotonic, objective 
function which may be used in place of the Gram ma- 
trix purity. After the purifications have been chosen, the 
benchmarking proceeds as in the pure state case (see Sec. 
II A) , requiring only the associated Gram matrix. 



D. Benchmarking with phase-symmetric ensembles 

One useful application for the above results is in testing 
devices which are phase covariant. A device (represented 
by A) is phase covariant when it commutes with unitary 
rotations of phase space, i.e., 

A[EVf/t] = U 4 A[a}Ul V a. (27) 

The rotation unitary is given by U$ — exp(—i<ph), with 
h the standard number operator, so t/J = U-^. We note 
that any channel can be made phase covariant by phase- 
randomization, 

Ac[ • ] = ^JJu 4 A[U-^ ■ )U^]Uld<j>. (28) 

One way to accomplish this phase randomization is to use 
a drifting optical phase [19], a situation which is com- 
mon in many continuous-variable setups. If the phase- 
randomized version of a channel passes a quantum bench- 
mark, then the channel itself must be in the quantum 
domain as well. On the other hand, since Ac involves a 
concatenation of channels, a phase- randomized channel 
may perform weaker against our benchmarks than the 
original channel. 

Phase-covariant channels offer a number of advantages. 
On the experimental side, they can be benchmarked us- 
ing only one physical test state. The effects of a phase- 
covariant channel on many other test states can be in- 
ferred by symmetry. On the theoretical side, phase co- 
variance can lead to great numerical simplifications. If 
the test states are unrelated, the number of free param- 
eters scales quadratically with the number of states M. 
When phase symmetry conditions are imposed, we will 
show that the number of free parameters can be made to 
scale only linearly with M. This gives us more compu- 
tational room to push up the number of states M in the 
test ensemble, leading to stronger benchmarks. 

To benchmark a phase-covariant channel, we consider 
test states from the rotationally-symmetric ensemble 

{pk\ Pk = Ul P0 Ul k }^\ (29) 

where 6 = |f (the angle 9 will always be defined this 
way). For phase-covariant devices, we have 

A[p fc ] = U*k[ Po \ul k , (30) 

so we can infer the channel's effect on p k by suitable rota- 
tions of its action on one real test state pq. In this way we 
can generate multi-state benchmark data efficiently from 
any seed state. We will refer to the situation where both 
the test states are phase symmetric and the tested device 
is phase covariant as phase-symmetric benchmarking. 

We will now present two standard form results relevant 
to phase-symmetric benchmarking. Given a particular 
device (mapping from system A' A" to B), the output 
state p°ab is defined by 

p% := [id A ® A] |* ent >(* ent UM" • ( 31 ) 
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where \^ cnt ) aa' A" ^ s the effective entangled state in Eq. 
(7). We must show that this reduced state is entangled 
in order to certify the given device. Without loss of gen- 
erality, we can decompose this state as 

M-lM-l 

Pab=mY, E \W\a®Pm\ (32) 
fe=0 z=o 

where p™* are square matrices on the B subsystem (the 
diagonal blocks p°^ — p£ ut are the output states). In a 
phase-symmetric benchmarking situation, the following 
relation can be imposed: 

UePkiUl = Pk + i,i + i (mod M) V k, I. (33) 

States with this symmetry can be brought into a simple 
standard form: 

Theorem 2 (Standard form) Let 

^ M-l M-l 

Tab =mY1 E \ k )( l \A®r kl , (34) 

k=0 1=0 

be an arbitrary bipartite matrix. Let 9 = and let ojm = 
exp(i#) denote the primitive Mth root of unity. Assume 
that the following symmetry relation holds: 

UeTkiUl = Tk+i.i+i {mod M) V k, I. (35) 

Then tab is unitarily equivalent to a block diagonal ma- 
trix 

M-l 

V(t AB ) = E k (36) 

fe=0 

where 

M-l 

Ek = m E "mthUI (37) 

1=0 

Moreover, 

M-l 
k=0 

and tab is positive semidefinite if and only if all of the 
Ek are positive semidefinite. 

Proof. See Appendix B. 

In the next section, we will quantify the entanglement 
of the state p A u B by computing lower bounds on the neg- 
ativity [35-38] , an entanglement measure defined by 

M{t ab ) = ll^lli-^*) , (39) 

where Ta denotes partial transposition and 1 1 • 1 1 1 is the 
trace norm. For phase-symmetric benchmarking, we can 
greatly simplify this calculation. 



Theorem 3 (Trace norm in standard form) Let 

tab be as in Theorem 2. Then the trace norm of the 
partially transposed state rj^ reduces to the form 

M-l 

\\t1 A b\\i= E II^Hi. ( 4 °) 

fe=0 

where the {-Efc}^} 1 are formed by rearranging the matrix 
elements of {Ek} 1 ^^ f rom the standard form of tab- 
Specifically, the matrix elements in the Fock basis are de- 
termined via 

[E k ]jl = [Ej+l-k)jl {mod M). (41) 

Proof. See Appendix C. 

These two standard form results allow us reduce the 
numerical complexity of phase-symmetric benchmarking. 
Although the constraints remain unchanged (the Gram 
matrix pa and measurements on the seed state po), the 
benchmarking state is more compactly encoded by using 
the M square matrices E k instead of the M ( M+1 ^ inde- 
pendent blocks tm- 



III. DEMONSTRATION WITH REAL MEMORY 
DATA 

We now demonstrate the above theoretical results us- 
ing real data generated by a continuous-variable mem- 
ory. This memory is based on a three-level lambda gra- 
dient echo scheme involving hyperfinc spin states of warm 
rubidium atoms. Its fundamental operation has been 
outlined in detail elsewhere [25]. The data used in the 
present paper is based on the results reported in [24]. 
This memory is not expected to have any sensitivity to 
phase. Indeed, the equations governing the memory oper- 
ation are phase- independent [25] . As well, there is no ex- 
perimental indication of phase sensitivity in the memory, 
i.e., the observed loss and excess noise are uniform with 
respect to the angle in phase space. Furthermore, there 
is no cross-talk between consecutive runs of the experi- 
ment. A data set contains measurements from 100,000 
pulses stored in, and then recalled from, the memory. 
Each pulse was separated in time by 0.1 ms which, given 
the decoherence rate of the atoms in the hot gas cell, 
means that each pulse was stored in a fresh ensemble of 
atoms. 

The experimental phase reference was allowed to drift 
freely between runs. In principle, this drift could be used 
to phase-randomize the device following the approach 
of [19], rendering the channel phase covariant. For the 
experimental data used in the present analysis, the ob- 
served drift was not fast enough to ensure the uniformly 
random distribution of phase angles required in Eq. (28). 
In future experiments, uniformity of the phase drift could 
be imposed by modifying the repetition rate of the ex- 
periment (see Sec. Ill A for more specific details). 



However, given the very strong evidence that the mem- 
ory is phase-insensitive and that consecutive runs are un- 
corrected, it is reasonable to expect that the memory 
operation is intrinsically phase-covariant. Our goal here 
is to demonstrate the strength of our benchmarking ap- 
proach under a variety of scenarios, by using data from 
a real memory. We therefore assume that running the 
experiment with a completely uniform phase drift would 
lead to the same data as obtained at the current rep- 
etition rate. In other words, the obtained experimental 
data is assumed to come from a phase-covariant memory. 
The phase-covariance assumption, which is independent 
from the benchmarking theory, allows us to easily gener- 
ate a number of different test ensembles using data from 
a single seed state, i.e., we have many non-orthogonal 
test states {pk\pk — UqPqUq To fully certify the 

memory with no assumptions (i.e., where we treat the 
memory as a black box) would require a more uniform 
phase drift than we obtained. Nevertheless, the current 
scenario allows for a useful demonstratation of our bench- 
marking tools. 



A. Experimental data acquisition 

During an experimental run, light is generated from 
a continuous-wave master laser and split into a coupling 
field (for activating the memory), a strong continuous- 
wave local oscillator, a weak signal pulse, and a memory- 
reference pulse. The memory-reference pulse is detuned 
from the active frequencies of the memory (while main- 
taining a well-defined phase relation with the signal 
pulse) , so that it passes through the memory without ab- 
sorption. The signal pulse is sent to the memory shortly 
after the reference pulse such that the phase drift be- 
tween the two pulses can be neglected. Thus, we can 
obtain phase information for the input signal by deter- 
mining the phase of the memory-reference pulse. 

Both the output signal and memory-reference pulses 
are interfered with the local oscillator (which has not 
passed through the memory) in a homodyne measure- 
ment scheme. Data was collected using a 12-bit resolu- 
tion data acquisition system. The homodyne measure- 
ment angle is varied between runs by allowing the signal 
and local oscillator to freely drift. This phase angle is 
identified by fitting a sine wave to the memory-reference 
pulse. A quadrature value was assigned to each signal 
pulse after interference with the local oscillator by inte- 
grating the amplitude of the pulse over its duration. Note 
that any inherent frequency shift due to the memory was 
cancelled coherently by adjusting the external magnetic 
field applied to the memory before data collection. This 
process generates quadrature measurement data for the 
output state pg ut . Similar data is obtained for the in- 
put state po by measuring it after it passes through an 
inactive memory [40]. 

Over repeated runs, many data points are collected to 
form a raw data set covering many phase angles. To en- 



Operator 


First moment (•} 


Second moment (■) 


X 


0.01 ± 0.03 


0.57 ± 0.04 


P 


-0.95 ± 0.03 


1.41 ± 0.09 



TABLE I: Experimentally-determined first and second mo- 
ments for the conjugate quadratures x and p. The vacuum 
state would have variances Var(a;) = Var(p) = | in these 
units. 



sure complete randomisation of the input pulse, the drift 
resulting from thermal fluctuations needs to be much 
faster than the repetition rate of the experiment. In our 
experimental setup, the pulse's thermal drift occurred 
at a rate on the order of Hz, while the experimental 
data was taken at a kHz rate. Running the experiment 
at a slower repetition rate would allow the phase refer- 
ence enough time to fully randomize, making the phase- 
randomization procedure of [19] possible and guarantee- 
ing phase-covariance. Of course, this would also affect 
the time needed to acquire the experimental data. 

Tomographic reconstruction was performed using an 
iterative maximum-likelihood algorithm [41] using data 
from 100,000 pulses, giving density matrices for both 
the input and output states in the Fock basis. A con- 
stant memory-induced phase offset [42] between the in- 
put and output signals was observed, and cancelled digi- 
tally during the tomography step. The tomographically- 
reconstructed density matrices were numerically trun- 
cated after the first 30 Fock states (|0) , . . . , |29)). This 
cutoff is supported by the fact that both the input state 
and output state are essentially confined to the first 10 
Fock levels, for which the cutoff error in either state 
is w 10 -5 . From the tomography, the input state was 
found to have mean photon number (nj n ) = 0.67 [24]. 
This value is consistent with the value obtained using 
direct quadrature measurement [43]. We do not ex- 
plicitly categorize this state, but it is qualitatively a 
'coherent state with added noise'. This extra noise is 
primarily due to amplitude fluctuations and changes in 
mode-matching occuring during the testing process. The 
tomographically-reconstructed output state had a similar 
form, with mean photon number (n ou t) = 0.57. 

As stated in Sec. II A, output state tomography is not 
required in quantum benchmarking; measurement of two 
conjugate quadratures may be sufficient. To compare the 
two approaches, we also consider benchmarks using mea- 
surements from only x and p. To get the required expec- 
tation values, we collected the original quadrature/phase 
data into bins containing 500 data points and restricted 
our attention to two bins corresponding to phase angles 
of and 90 degrees. Taking the mean and standard de- 
viation of the data within these bins, we get estimates of 
(x), (p), Var(x), and Var(p), where Var(z) = (i 2 ) — (z) 2 
for an operator z £ {x,p}. Error bars, based on the 
finite sample size of the bins, are calculated using stan- 
dard error propagation techniques. The final quadrature 
moment values are listed in Table I. 



9 



B. Finding the Gram matrix 

The test ensemble is generated by rotations of the seed 
state po, as in Eq. (29), forming a ring in phase space. 
Our theoretical description of the seed state comes from 
the tomographic reconstruction with the finite cutoff in 
Fock space. We make no further assumptions about the 
seed state other than that the tomographic description is 
accurate. Purifications, and the associated Gram matrix, 
will be based on this numerical description. Due to the 
mixed nature of this state, previous entanglement-based 
benchmarking schemes would be limited to only two test 
states. We will see later that two test states are not 
sufficient in this case to reliably (i.e., within experimen- 
tal error) give non-zero output entanglement, but adding 
more test states enables a successful benchmarking. 

In order to find candidate Gram matrices for the pu- 
rifications, we convert the optimization to a semidefinite 
program. To this end, we consider the bipartite matrix 

pJ A , =Tr A „ |* ent >(* eIlt U<A'< 

M-1M-1 



(42) 



fe=0 1=0 



The diagonal blocks of p l AA , are the test states, i.e., 
PTk = Pk i an d the off-diagonal blocks p' fe " consist of free 
parameters. The off-diagonal blocks are linked to the 
purifications through the relation 



Tr(p%)=M(k\p A \l) = (T l \T k ) 



(43) 



We have traced out the purifying system A" since no 
measurements are performed on that system, but we have 
kept subsystem A' in order to later ascribe some quantity 
of entanglement to the input state. 

As stated in Sec. II C, instead of using the Gram ma- 
trix purity V, we optimized the slightly different and nu- 
merically more simple objective function h from Eq. (26). 
Although it might be more elegant to use the Gram ma- 
trix purity V , nothing in the benchmarking procedure de- 
pends on that function. In fact, the Gram matrix found 
by optimizing h over the given experimental data works 
quite well for our purposes. To see this, we can consider 
simple lower and upper bounds on the maximum value 
of V compatible with the given constraints. 

For one, the purity of the Gram matrix that maximizes 
h (denoted arg max h) is less than the maximum possible 
Gram matrix purity over all purifications, i.e., 



V I arg max h < max7>. 



(44) 



On the other hand, we can easily find an upper bound 
on V by considering the pairwise fidelities: 



M-XM-X 

maxP < —j F iPk,Pi)- 



(45) 



fc=0 1=0 



In Fig. 1, we plot the values of these lower and upper 
bounds for different numbers of rotationally symmetric 



^ 1 



0.4 



0.35 



0.3 



.5 J3 0.25 



O £ 



0.2 



-T- Upper bound 
--»- Lower bound 



0.4 
0.39 
0.38 
0.37 
0.36 



0.15 

2 3 4 5 6 7 8 

Number of test states M 

FIG. 1: Upper and lower bounds for the maximum purity V of 
experimentally-compatible Gram matrices, corresponding to 
rotationally symmetric test ensembles of different cardinality. 
We have removed the factor , which is common to both the 
upper and lower bounds, and independent of the purifications. 
The difference between the two bounds is on the order of 
10 -3 or less (inset), meaning that the obtained Gram matrix 
has purity which is near-optimal (represented by the lower 
curve). The estimated numerical precision of the optimization 
is « 10" 3 . 



test states (omitting the constant factors jjj)- The com- 
puted bounds are very close to each other, differing only 
in the third decimal place, which is the same order as 
the estimated numerical precision of our optimization. 
We conclude that using the modified objective function 
in Eq. (26) is suitable for the situation we are consider- 
ing, and that the obtained Gram matrix purity is close 
to what would be obtained by maximizing V itself. 



C. Computing the entanglement 

Having used the heuristic function h to find a feasible 
P'aa' which corresponds to a near-optimal Gram matrix, 
we now consider the entanglement of the corresponding 
output state p° A ^ . We use the negativity as our entangle- 
ment measure, which can be computed as a semidefinite 
program, namely 



N{t ab ) = min Tr(r_) 

T-£T 



(46) 



where the constraint set T is defined by 

T := {r_ |r_ > and t a % + r_ > 0}. (47) 

For optimizations where t A b is not fully known (as is 
the situation here) , constraints from experimental obser- 
vations are also included in this set. For states written 
in standard form, the negativity can be computed on the 
level of the individual matrices E k . 

Next, we used the entanglement quantification tools 
from [21] to compute the lowest value of the negativ- 
ity that could be compatible with the given data. In 
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this approach, one rigorously finds lower bounds on the 
negativity of the output state p°^g by truncating the B 
subsystem of the density matrix at some finite Fock level 
N , which is large compared to the observed mean pho- 
ton number n. For this test, we imposed the cutoff at 
Fock level N = 15, making the truncated matrices 16- 
dimensional. 

The quantification procedure was applied to four mea- 
surement scenarios: 

(i) Tomography: the output density matrix is fully 
known numerically (up to the truncation level N); 

(ii) Quadratures from tomography: only mea- 
surements on conjugate quadratures x and p for 
the output state are used, and these values are 
obtained from the tomographically-reconstructed 
density matrix; 

(iii) Quadratures from data: same as (ii), but the 
quadrature values are obtained directly from the 
homodyne data; 

(iv) Quadratures with error bars: the same quadra- 
ture data is used as in (iii), but estimates of data 
error bars are included in the computation (o — 3cr 
levels). 

In cases (ii)-(iv), the choice of TV affects the strength of 
the optimization procedure itself [21], not just the den- 
sity matrix. Truncating at a higher level would lead to 
better bounds on the entanglement. In all cases, we use 
the standard form results of Theorems 2 and 3. Finally, 
although [21] does not explicitly give a procedure for 
including error bars in the computation, it is straight- 
forward to relax the framework and include them. We 
performed the optimization in Matlab, using the fron- 
tcnd YALMIP [44] and the solver SDPT3 [45]. The nu- 
merical precision of the optimization is estimated to be 
rj 10~ 3 . The results of these computations for different 
sized rotationally-symmetric ensembles is shown in Figs. 
2 and 3. 

There are a number of things to point out about the 
results in Figs. 2 and 3. Primarily, we see that our 
benchmarking tools were successful in finding non-zero 
output entanglement in this phase-symmetric situation. 
Not only is the output state p°Jjj entangled for all M 
when we use the tomographically reconstructed density 
matrix /?Q Ut , but we can also see non-zero entanglement 
using the restricted set of measurements, including error 
bars. Furthermore, there is an evident advantage in using 
more than two test states, since it leads to higher entan- 
glement in all cases. Indeed, it is not until M > 2 that 
the 3er level conclusively leads to non-zero entanglement 
in the output state. 

Another interesting result is that scenarios (i) and 
(ii) give very similar results for the entanglement when 
M > 3. This is likely due to the fact that the output state 
has low levels of excess noise. The quadratures would 
already reveal that such a state is close to a minimum 



--«- Input state 

Output state (tomo.) 
- -*- Output state (quad, from tomo.) 




2 3 4 5 6 7 

Number of test states M 



FIG. 2: Benchmarking results based on rotationally- 
symmetric test ensembles using tomographically- 
reconstructed density matrices. The top curve is the 
negativity of the input state P^a' found by optimizing the 
heuristic function h. The other curves are lower bounds 
to the negativity of p A n B based on (in descending order): 
(i) tomographic reconstruction and (ii) quadrature values 
obtained from the tomographic reconstruction. A bench- 
marking procedure is successful when the output state has 
negativity larger than zero. 



Input state 

Output state (quad, from data) 




2 3 4 5 6 7 

Number of test states M 



FIG. 3: Benchmarking results based on rotationally- 
symmetric test ensembles using direct quadrature measure- 
ment data. The top curve is the same as Fig. 2. The other 
curves are lower bounds to the negativity of p°ab based on (in 
descending order): (iii) quadrature values obtained directly 
from measurement data and (iv) quadrature data with er- 
ror bars included in the optimization (la, 2a, and 3a levels). 
For more than 3 test states, there is non-zero entanglement 
remaining even when including error estimates. 



uncertainty state. Additional measurement information 
therefore does not contribute significantly to our knowl- 
edge of the state. For other types of states, quadrature 
measurements and tomography should not be expected to 
yield similar results (though there may be other choices 
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of measurement operators which work well). Note that 
we do not compare situations (i) and (iii) because of the 
extra phase shift included in the tomographic reconstruc- 
tion outlined in Sec. Ill A. 

We also point out the advantage offered by phase- 
symmetric benchmarking and the standard form. Pre- 
viously, entanglement quantification results were demon- 
strated for test ensembles consisting of only two or three 
states [20, 21]. Although the methods of [21] apply to 
any number of (pure) test states, the computational re- 
sources required for more than three test states were pro- 
hibitive. Here, using the standard form, we have no prob- 
lem pushing up to much higher numbers of states; Figs. 
1-3 contain results for up to 8 test states. In fact, with 
a desktop PC, we were able to calculate benchmarks for 
up to 10 states, but these results are not included in 
the figures because they do not significantly change af- 
ter M = 8. Finally, we note that all the results for the 
output state entanglement are lower bounds. If we were 
to use a real entangled state instead of a virtual one, 
we would have access to more information, and the out- 
put entanglement would be higher. Even still, for bench- 
marking schemes based on a virtual entangled state, the 
entanglement quantification procedure yields useful and 
illuminating results. 



IV. CONCLUSION 

We have extended the methodology for entanglement- 
based quantum benchmarking to arbitrary finite ensem- 
bles of test states. This extension considers purifications 
of the test states, and we have outlined how to find strong 
choices for these purifications, leading to the best chance 
of certifying a device in the quantum domain. We also 
gave equations for reducing the test states to a standard 
form when certain phase symmetry conditions hold, al- 
lowing us to generate benchmarks from any single seed 
state. These theoretical tools were demonstrated using 
data from a real-world quantum memory and to quan- 
titatively explore different benchmarking scenarios. Al- 
though we mainly considered test states which were of a 
continuous-variable nature, all of our results, except the 
standard form, apply to discrete states as well. Indeed, 
the same framework can be used to calculate fidelity- 
based benchmarks, which may be better suited to low- 
dimensional states [46]. Since our benchmarking scheme 
works for arbitrary test states, it incorporates previous 
entanglement-based benchmarking tools [15-21] into a 
unified framework. Together, these tools allow us to per- 
form efficient benchmarking tests on quantum devices us- 
ing realistic experimental resources. 
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Appendix A: Convexity proofs 

In this appendix, we prove the following two state- 
ments about convexity: i) the objective function in Eq. 
(19) is a convex function, and ii) the set of Gram matrices 
which are consistent with purifications of the test states 
is a convex set. First, for positive semidefinite operators, 
the purity is equivalent to the squared Hilbert-Schmidt 
norm. By definition, all norms are convex, and taking 
the square preserves this convexity. Hence, our objective 
function is convex. 

To prove the second statement, let {|r£) ,„} fc= ^ 1 and 
{Irfc)^/^"}/^ 1 be- two sets of purifications of the test 
states {/5/cjfclo 1 (without loss of generality, we can con- 
sider the purifying system to be the same), with Gram 
matrices Go and G\, respectively. Fix some p € [0,1] 
and take the convex combination of Gram matrices G = 
pGo + (1 — p)G\. We need to show that there is some 
compatible set of purifications leading to the Gram ma- 
trix G. 

To do this, we define the following states: 

\Xk) A >A»A»> ^ VP" K) A , A „ |0> a ,„ + V^P K) A >A» |1>A'« > 

(Al) 

where |0)^,„ and |1)^,// are orthonormal states on some 
additional purifying system A'". Tracing out all systems 
except A' , we find 

Tr A " A 'n \Xk){Xk\ A >A"A"< = P Tv A" | r fe)( r feL'A" 

+ (l-p)TT A "\rl)(Tl\ A , At , 

= Pk- (A2) 

Thus, the \xk) a'A"A"' are purifications of the pk- The 
elements of the corresponding Gram matrix are 

(xk\x l ) = p{n\r { !) + (i-p) (rl\rj) 

= p[G ]fci + (l-p)[Gi] w , (A3) 

which are exactly the elements of the convex combination 
G. Therefore, the set of Gram matrices which come from 
purifications of the test states is a convex set. 



Appendix B: Proof of standard form 

In this appendix, we give the proof of the standard 
form in Theorem 2. If Eq. (35) holds, then t A b can be 
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transformed by the unitary matrix 

M-l 



R=®u$ 



(Bl) 



k=0 



to the form 
C(t ab ) = R)t A bR 



1 

M 



Too Wqi W 2 

Wo, M-l TOO Wqi 
Wq.M-2 Wo, M-l T 00 



w, 



01 



w 



02 



W 



O.i 



W), M-l 
Wo,M-2 
W M-3 



TOO 



(B2) 



with Wy := TijUg. The matrix C(tab) has a 6/ocfc cir- 
culant structure, meaning that each row of blocks is the 
same as the previous row, but shifted by one to the right. 
Now, let u>m — exp(i9) be the primitive Mth root of 
unity and let 



i , 

7m M 



(B3) 



be the (unitary) discrete Fourier transform matrix. Us- 
ing a standard theorem on block circulant matrices [47], 
C(tab) can be block diagonalized by the matrix Fm<8>1b- 
Explicitly, 



V(t ab ) =(Fl ® 1 b )C(t ab )(F m ® Is) 

M-l 

= E k , (B4) 



k=0 



with 



M-l 

E * = M E W ™ 



(B5) 



(=0 



Since the partial trace is unaffected by local unitaries on 
the same subsystem, we have 



M-l /M-l 

]T E k =Tr A l fe X*L 
fc=0 V fe=0 / 

=Tr^C(rAB). 
Comparing this with Eq. (B2), we conclude 



(B6) 



Appendix C: Proof of trace norm in standard form 

In this appendix, we prove Theorem 3, which gives a 
formula for the negativity involving the standard form. 
If tab satisfies the symmetry condition in Eq. (35), 
then so will the partial transpose rj^. Therefore, rj^ 
is unitarily equivalent to some block diagonal matrix 
B( t ab) = ®k=o~ -^h- The trace norm is a unitarily in- 
variant norm, so we must have 

M-l 



l^lll 



0^ 

fc=0 



M-l 



k 1, 



(CI) 



fc=0 



which proves Eq. (40). 

In order to determine the matrices E k in Eq. (41), 
we first find an explicit expression relating the standard 
form matrices E k and the matrix tab ■ Define the matrix 
elements of a bipartite matrix H by 



ij,kl 



(j\ B H\k) A <E>\l) 1 



(C2) 



For subsystem B, we work in the Fock basis, where the 
rotation operator Ug is diagonal, with elements 



[U e ]ji = uJ s ji- 



(C3) 



The unitary matrix R in Eq. (Bl), formed by taking 
powers of Ug, is also diagonal in this basis, with the fol- 
lowing elements: 



[R]ij,kl 



■>M 3 &ik8jl- 



(C4) 



Therefore, the Fock basis elements of the block circulant 
matrix C(tab), defined in Eq. (B2), are given by 



[C(TAB)]ij,kl = W M [TAB]ij,kl 



(C5) 



To get the standard form, we need to perform a Fourier 
transform on the block circulant matrix C(tab) as in Eq. 
(B4) . The elements of the discrete Fourier transform Fm 
are given in Eq. (B3). After substitution, we arrive at 
the matrix elements of the standard form: 



P{rA B )} 



ij,kl 



^ M-l M-l 

mE £< 

m— n— 
A/-1M-1 
1 \ " \ " m(j-i)+n{k-l) 
Z M 1^ 1^ w m 

m — n—O 



my.nl 



T AB\mj 1 nl' 

(C6) 



M-l 



E^ = 



TOO, 



fc=0 



(B7) 



which proves Eq. (38). Finally, the positive semidefinitc 
condition follows directly from the unitary equivalence of 
tab and V(t A b)- 

□ 



Since V{tab) is block diagonal, V(tab) 
have 



/tnM-1 
U7fc=0 



Ek, we 



[E k }ji :=[D{t A b)] 



kj,kl 



^ M-l M-l 

m E E< 

m— n— 



m(j — k)-\-n(k—l) 
J M 



T AB\mj,nl- 



(C7) 
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This is the final formula linking the elements of tab and 
the matrices {Ek}^ 1 ^ 1 ■ For completeness, the inverse 
formula is given by 



Hence, the the matrix E indexed by (J + I — k) has ele- 
ments 



[TAB\ij,H - jj u m 



M-l 

1 i rn(i — k)+k-l — i-j 

m=0 



E 



m\jl- 



(C8) 



When is used in this formula instead of tab, we 
find the relation 



M-l M-l 

m(j-k)+n(k-l), i /p Q \ 
M ^ ^ M [TAB\nj,rnl- (C9) 

m— n— 



_^ M—l M— ] 



M-l M-l 



m(k — l)-\-n(j — k) 



M ^ ^ ~ M 

m—0 7i—0 



(CIO) 

which is the same as Eq. (C7), after interchanging the 
summation indices m and n. Therefore, [Ej+i—k]ji — 
[Ek]ji (in the Fock basis). 
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